Salivary gland carcinosarcoma ex pleomorphic adenoma; a case report with molecular characterization.

True malignant mixed tumors, also known as salivary gland carcinosarcoma (SCS), are uncommon yet highly aggressive lesions associated with a poor prognosis. These tumors exhibit a distinctive biphasic structure characterized by both epithelial and mesenchymal components. Recent research has shown that the majority of SCS cases stem from pre-existing pleomorphic adenomas (PAs), suggesting a stepwise developmental pattern. In this report, we present a case of a 73-year-old female with SCS and describe the clinical, radiographic, and pathologic observations. Notably, the SCS was associated with a residual PA. The SCS displayed a CTNNB1::PLAG1 gene rearrangement, providing a molecular basis for its origin from the PA. Further DNA genomic analysis exposed mutations in BAP1, PER1, and LRPB1. Our findings provide support to the theory that SCS emerges from a pre-existing PA while highlighting the multiple genetic changes that could contribute to malignant transformation.

Increase in Full-Length Dystrophin by Exon Skipping in Duchenne Muscular Dystrophy Patients with Single Exon Duplications: An Open-label Study.

Single exon duplications account for disease in a minority of Duchenne muscular dystrophy patients. Exon skipping in these patients has the potential to be highly therapeutic through restoration of full-length dystrophin expression. We conducted a 48-week open label study of casimersen and golodirsen in 3 subjects with an exon 45 or 53 duplication. Two subjects (aged 18 and 23 years) were non-ambulatory at baseline. Upper limb, pulmonary, and cardiac function appeared stable in the 2 subjects in whom they could be evaluated. Dystrophin expression increased from 0.94 % ±0.59% (mean±SD) of normal to 5.1% ±2.9% by western blot. Percent dystrophin positive fibers also rose from 14% ±17% at baseline to 50% ±42% . Our results provide initial evidence that the use of exon-skipping drugs may increase dystrophin levels in patients with single-exon duplications.

Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity.

BACKGROUND: Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms. RESULTS: In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms. CONCLUSION: This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

CRISPR-Cas9 homology-independent targeted integration of exons 1-19 restores full-length dystrophin in mice.

Duchenne muscular dystrophy is an X-linked disorder typically caused by out-of-frame mutations in the DMD gene. Most of these are deletions of one or more exons, which can theoretically be corrected through CRISPR-Cas9-mediated knockin. Homology-independent targeted integration is a mechanism for achieving such a knockin without reliance on homology-directed repair pathways, which are inactive in muscle. We designed a system based on insertion into intron 19 of a DNA fragment containing a pre-spliced mega-exon encoding DMD exons 1-19, along with the MHCK7 promoter, and delivered it via a pair of AAV9 vectors in mice carrying a Dmd exon 2 duplication. Maximal efficiency was achieved using a Cas9:donor adeno-associated virus (AAV) ratio of 1:5, with Cas9 under the control of the SPc5-12 promoter. This approach achieved editing of 1.4% of genomes in the heart, leading to 30% correction at the transcript level and restoration of 11% of normal dystrophin levels. Treatment efficacy was lower in skeletal muscles. Sequencing additionally revealed integration of fragmentary and recombined AAV genomes at the target site. These data provide proof of concept for a gene editing system that could restore full-length dystrophin in individuals carrying mutations upstream of intron 19, accounting for approximately 25% of Duchenne muscular dystrophy patients.

Multiomics in primary and metastatic breast tumors from the AURORA US network finds microenvironment and epigenetic drivers of metastasis.

The AURORA US Metastasis Project was established with the goal to identify molecular features associated with metastasis. We assayed 55 females with metastatic breast cancer (51 primary cancers and 102 metastases) by RNA sequencing, tumor/germline DNA exome and low-pass whole-genome sequencing and global DNA methylation microarrays. Expression subtype changes were observed in ~30% of samples and were coincident with DNA clonality shifts, especially involving HER2. Downregulation of estrogen receptor (ER)-mediated cell-cell adhesion genes through DNA methylation mechanisms was observed in metastases. Microenvironment differences varied according to tumor subtype; the ER+/luminal subtype had lower fibroblast and endothelial content, while triple-negative breast cancer/basal metastases showed a decrease in B and T cells. In 17% of metastases, DNA hypermethylation and/or focal deletions were identified near HLA-A and were associated with reduced expression and lower immune cell infiltrates, especially in brain and liver metastases. These findings could have implications for treating individuals with metastatic breast cancer with immune- and HER2-targeting therapies.

Pacific Biosciences Fusion and Long Isoform Pipeline for Cancer Transcriptome-Based Resolution of Isoform Complexity.

Genomic profiling using short-read sequencing has utility in detecting disease-associated variation in both DNA and RNA. However, given the frequent occurrence of structural variation in cancer, molecular profiling using long-read sequencing improves the resolution of such events. For example, the Pacific Biosciences long-read RNA-sequencing (Iso-Seq) transcriptome protocol provides full-length isoform characterization, discernment of allelic phasing, and isoform discovery, and identifies expressed fusion partners. The Pacific Biosciences Fusion and Long Isoform Pipeline (PB_FLIP) incorporates a suite of RNA-sequencing software analysis tools and scripts to identify expressed fusion partners and isoforms. In addition, sequencing of a commercial reference (Spike-In RNA Variants) with known isoform complexity was performed and demonstrated high recall of the Iso-Seq and PB_FLIP workflow to benchmark our protocol and analysis performance. This study describes the utility of Iso-Seq and PB_FLIP analysis in improving deconvolution of complex structural variants and isoform detection within an institutional pediatric and adolescent/young adult translational cancer research cohort. The exemplar case studies demonstrate that Iso-Seq and PB_FLIP discover novel expressed fusion partners, resolve complex intragenic alterations, and discriminate between allele-specific expression profiles.

Molecular Heterogeneity in Pediatric Malignant Rhabdoid Tumors in Patients With Multi-Organ Involvement.

Rhabdoid tumors (RTs) of the brain (atypical teratoid/rhabdoid tumor; AT/RT) and extracranial sites (most often the kidney; RTK) are malignant tumors predominantly occurring in children, frequently those with SMARCB1 germline alterations. Here we present data from seven RTs from three pediatric patients who all had multi-organ involvement. The tumors were analyzed using a multimodal molecular approach, which included exome sequencing of tumor and germline comparator and RNA sequencing and DNA array-based methylation profiling of tumors. SMARCB1 germline alterations were identified in all patients and in all tumors. We observed a second hit in SMARCB1 via chr22 loss of heterozygosity. By methylation profiling, all tumors were classified as rhabdoid tumors with a corresponding subclassification within the MYC, TYR, or SHH AT/RT subgroups. Using RNA-seq gene expression clustering, we recapitulated the classification of known AT/RT subgroups. Synchronous brain and kidney tumors from the same patient showed different patterns of either copy number variants, single-nucleotide variants, and/or genome-wide DNA methylation, suggestive of non-clonal origin. Furthermore, we demonstrated that a lung and abdominal metastasis from two patients shared overlapping molecular features with the patient's primary kidney tumor, indicating the likely origin of the metastasis. In addition to the SMARCB1 events, we identified other whole-chromosome events and single-nucleotide variants in tumors, but none were found to be prognostic, diagnostic, or offer therapeutic potential for rhabdoid tumors. While our findings are of biological interest, there may also be clinical value in comprehensive molecular profiling in patients with multiple rhabdoid tumors, particularly given the potential prognostic and therapeutic implications for different rhabdoid tumor subgroups demonstrated in recent clinical trials and other large cohort studies.

High early death rates, treatment resistance, and short survival of Black adolescents and young adults with AML.

Survival of patients with acute myeloid leukemia (AML) is inversely associated with age, but the impact of race on outcomes of adolescent and young adult (AYA; range, 18-39 years) patients is unknown. We compared survival of 89 non-Hispanic Black and 566 non-Hispanic White AYA patients with AML treated on frontline Cancer and Leukemia Group B/Alliance for Clinical Trials in Oncology protocols. Samples of 327 patients (50 Black and 277 White) were analyzed via targeted sequencing. Integrated genomic profiling was performed on select longitudinal samples. Black patients had worse outcomes, especially those aged 18 to 29 years, who had a higher early death rate (16% vs 3%; P=.002), lower complete remission rate (66% vs 83%; P=.01), and decreased overall survival (OS; 5-year rates: 22% vs 51%; P<.001) compared with White patients. Survival disparities persisted across cytogenetic groups: Black patients aged 18 to 29 years with non-core-binding factor (CBF)-AML had worse OS than White patients (5-year rates: 12% vs 44%; P<.001), including patients with cytogenetically normal AML (13% vs 50%; P<.003). Genetic features differed, including lower frequencies of normal karyotypes and NPM1 and biallelic CEBPA mutations, and higher frequencies of CBF rearrangements and ASXL1, BCOR, and KRAS mutations in Black patients. Integrated genomic analysis identified both known and novel somatic variants, and relative clonal stability at relapse. Reduced response rates to induction chemotherapy and leukemic clone persistence suggest a need for different treatment intensities and/or modalities in Black AYA patients with AML. Higher early death rates suggest a delay in diagnosis and treatment, calling for systematic changes to patient care.

Exome sequencing in multiplex families with left-sided cardiac defects has high yield for disease gene discovery.

Congenital heart disease (CHD) is a common group of birth defects with a strong genetic contribution to their etiology, but historically the diagnostic yield from exome studies of isolated CHD has been low. Pleiotropy, variable expressivity, and the difficulty of accurately phenotyping newborns contribute to this problem. We hypothesized that performing exome sequencing on selected individuals in families with multiple members affected by left-sided CHD, then filtering variants by population frequency, in silico predictive algorithms, and phenotypic annotations from publicly available databases would increase this yield and generate a list of candidate disease-causing variants that would show a high validation rate. In eight of the nineteen families in our study (42%), we established a well-known gene/phenotype link for a candidate variant or performed confirmation of a candidate variant's effect on protein function, including variants in genes not previously described or firmly established as disease genes in the body of CHD literature: BMP10, CASZ1, ROCK1 and SMYD1. Two plausible variants in different genes were found to segregate in the same family in two instances suggesting oligogenic inheritance. These results highlight the need for functional validation and demonstrate that in the era of next-generation sequencing, multiplex families with isolated CHD can still bring high yield to the discovery of novel disease genes.

Detection of brain somatic variation in epilepsy-associated developmental lesions.

OBJECTIVE: Epilepsy-associated developmental lesions, including malformations of cortical development and low-grade developmental tumors, represent a major cause of drug-resistant seizures requiring surgical intervention in children. Brain-restricted somatic mosaicism has been implicated in the genetic etiology of these lesions; however, many contributory genes remain unidentified. METHODS: We enrolled 50 children who were undergoing epilepsy surgery into a translational research study. Resected tissue was divided for clinical neuropathologic evaluation and genomic analysis. We performed exome and RNA sequencing to identify somatic variation and we confirmed our findings using high-depth targeted DNA sequencing. RESULTS: We uncovered candidate disease-causing somatic variation affecting 28 patients (56%), as well as candidate germline variants affecting 4 patients (8%). In agreement with previous studies, we identified somatic variation affecting solute carrier family 35 member A2 (SLC35A2) and mechanistic target of rapamycin kinase (MTOR) pathway genes in patients with focal cortical dysplasia. Somatic gains of chromosome 1q were detected in 30% (3 of 10) of patients with Type I focal cortical dysplasia (FCD)s. Somatic variation in mitogen-activated protein kinase (MAPK) pathway genes (i.e., fibroblast growth factor receptor 1 [FGFR1], FGFR2, B-raf proto-oncogene, serine/threonine kinase [BRAF], and KRAS proto-oncogene, GTPase [KRAS]) was associated with low-grade epilepsy-associated developmental tumors. RNA sequencing enabled the detection of somatic structural variation that would have otherwise been missed, and which accounted for more than one-half of epilepsy-associated tumor diagnoses. Sampling across multiple anatomic regions revealed that somatic variant allele fractions vary widely within epileptogenic tissue. Finally, we identified putative disease-causing variants in genes not yet associated with focal cortical dysplasia. SIGNIFICANCE: These results further elucidate the genetic basis of structural brain abnormalities leading to focal epilepsy in children and point to new candidate disease genes.

Discovery of clinically relevant fusions in pediatric cancer.

BACKGROUND: Pediatric cancers typically have a distinct genomic landscape when compared to adult cancers and frequently carry somatic gene fusion events that alter gene expression and drive tumorigenesis. Sensitive and specific detection of gene fusions through the analysis of next-generation-based RNA sequencing (RNA-Seq) data is computationally challenging and may be confounded by low tumor cellularity or underlying genomic complexity. Furthermore, numerous computational tools are available to identify fusions from supporting RNA-Seq reads, yet each algorithm demonstrates unique variability in sensitivity and precision, and no clearly superior approach currently exists. To overcome these challenges, we have developed an ensemble fusion calling approach to increase the accuracy of identifying fusions. RESULTS: Our Ensemble Fusion (EnFusion) approach utilizes seven fusion calling algorithms: Arriba, CICERO, FusionMap, FusionCatcher, JAFFA, MapSplice, and STAR-Fusion, which are packaged as a fully automated pipeline using Docker and Amazon Web Services (AWS) serverless technology. This method uses paired end RNA-Seq sequence reads as input, and the output from each algorithm is examined to identify fusions detected by a consensus of at least three algorithms. These consensus fusion results are filtered by comparison to an internal database to remove likely artifactual fusions occurring at high frequencies in our internal cohort, while a "known fusion list" prevents failure to report known pathogenic events. We have employed the EnFusion pipeline on RNA-Seq data from 229 patients with pediatric cancer or blood disorders studied under an IRB-approved protocol. The samples consist of 138 central nervous system tumors, 73 solid tumors, and 18 hematologic malignancies or disorders. The combination of an ensemble fusion-calling pipeline and a knowledge-based filtering strategy identified 67 clinically relevant fusions among our cohort (diagnostic yield of 29.3%), including RBPMS-MET, BCAN-NTRK1, and TRIM22-BRAF fusions. Following clinical confirmation and reporting in the patient's medical record, both known and novel fusions provided medically meaningful information. CONCLUSIONS: The EnFusion pipeline offers a streamlined approach to discover fusions in cancer, at higher levels of sensitivity and accuracy than single algorithm methods. Furthermore, this method accurately identifies driver fusions in pediatric cancer, providing clinical impact by contributing evidence to diagnosis and, when appropriate, indicating targeted therapies.

Evidence of pioneer factor activity of an oncogenic fusion transcription factor.

Recent characterizations of pioneer transcription factors provide insights into their structures and patterns of chromatin recognition associated with their roles in cell fate commitment and transformation. Intersecting with these basic science concepts, identification of pioneer factors (PFs) fused together as driver translocations in childhood cancers raises questions of whether these fusions retain the fundamental ability to invade repressed chromatin, consistent with their monomeric PF constituents. This study defines the cellular and chromatin localization of PAX3-FOXO1, an oncogenic driver of childhood rhabdomyosarcoma (RMS), derived from a fusion of PFs. To quantitatively define its chromatin-targeting functions and capacity to drive epigenetic reprogramming, we developed a ChIP-seq workflow with per-cell normalization (pc-ChIP-seq). Our quantitative localization studies address structural variation in RMS genomes and reveal insights into inactive chromatin localization of PAX3-FOXO1. Taken together, our studies are consistent with pioneer function for a driver oncoprotein in RMS, with repressed chromatin binding and nucleosome-motif targeting.

Endogenous retrovirus envelope as a tumor-associated immunotherapeutic target in murine osteosarcoma.

Osteosarcoma remains one of the deadliest cancers in pediatrics and young adults. We administered two types of immunotherapies, oncolytic virotherapy and immune checkpoint inhibition, to two murine osteosarcoma models and observed divergent results. Mice bearing F420 showed no response, whereas those with K7M2 showed prolonged survival in response to combination therapy. K7M2 had higher expression of immune-related genes and higher baseline immune cell infiltrates, but there were no significant differences in tumor mutational burden or predicted MHC class I binding of nonsynonymous mutations. Instead, we found several mouse endogenous retrovirus sequences highly expressed in K7M2 compared with F420. T cell tetramer staining for one of them, gp70, was detected in mice with K7M2 but not F420, suggesting that endogenous retrovirus proteins are targets for the anti-tumor immune reaction. Given prior observations of endogenous retrovirus expression in human osteosarcomas, our findings may be translatable to human disease.

Gastroblastoma with a novel EWSR1-CTBP1 fusion presenting in adolescence.

Gastroblastomas are rare tumors with a biphasic epithelioid/spindle cell morphology that typically present in early adulthood and have recurrent MALAT1-GLI1 fusions. We describe an adolescent patient with Wiskott-Aldrich syndrome who presented with a large submucosal gastric tumor with biphasic morphology. Despite histologic features consistent with gastroblastoma, a MALAT1-GLI1 fusion was not found in this patient's tumor; instead, comprehensive molecular profiling identified a novel EWSR1-CTBP1 fusion and no other significant genetic alterations. The tumor also overexpressed NOTCH and FGFR by RNA profiling. The novel fusion and expression profile suggest a role for epithelial-mesenchymal transition in this tumor, with potential implications for the pathogenesis of biphasic gastric tumors such as gastroblastoma.

Synonymous variants that disrupt messenger RNA structure are significantly constrained in the human population.

BACKGROUND: The role of synonymous single-nucleotide variants in human health and disease is poorly understood, yet evidence suggests that this class of "silent" genetic variation plays multiple regulatory roles in both transcription and translation. One mechanism by which synonymous codons direct and modulate the translational process is through alteration of the elaborate structure formed by single-stranded mRNA molecules. While tools to computationally predict the effect of non-synonymous variants on protein structure are plentiful, analogous tools to systematically assess how synonymous variants might disrupt mRNA structure are lacking. RESULTS: We developed novel software using a parallel processing framework for large-scale generation of secondary RNA structures and folding statistics for the transcriptome of any species. Focusing our analysis on the human transcriptome, we calculated 5 billion RNA-folding statistics for 469 million single-nucleotide variants in 45,800 transcripts. By considering the impact of all possible synonymous variants globally, we discover that synonymous variants predicted to disrupt mRNA structure have significantly lower rates of incidence in the human population. CONCLUSIONS: These findings support the hypothesis that synonymous variants may play a role in genetic disorders due to their effects on mRNA structure. To evaluate the potential pathogenic impact of synonymous variants, we provide RNA stability, edge distance, and diversity metrics for every nucleotide in the human transcriptome and introduce a "Structural Predictivity Index" (SPI) to quantify structural constraint operating on any synonymous variant. Because no single RNA-folding metric can capture the diversity of mechanisms by which a variant could alter secondary mRNA structure, we generated a SUmmarized RNA Folding (SURF) metric to provide a single measurement to predict the impact of secondary structure altering variants in human genetic studies.

YAP1-FAM118B Fusion Defines a Rare Subset of Childhood and Young Adulthood Meningiomas.

Meningiomas are a central nervous system tumor primarily afflicting adults, with <1% of cases diagnosed during childhood or adolescence. Somatic variation in NF2 may be found in ∼50% of meningiomas, with other genetic drivers (eg, SMO, AKT1, TRAF7) contributing to NF2 wild-type tumors. NF2 is an upstream negative regulator of YAP signaling and loss of the NF2 protein product, Merlin, results in YAP overexpression and target gene transcription. This mechanism of dysregulation is described in NF2-driven meningiomas, but further work is necessary to understand the NF2-independent mechanism of tumorigenesis. Amid our institutional patient-centric comprehensive molecular profiling study, we identified an individual with meningioma harboring a YAP1-FAM118B fusion, previously reported only in supratentorial ependymoma. The tumor histopathology was remarkable, characterized by prominent islands of calcifying fibrous nodules within an overall collagen-rich matrix. To gain insight into this finding, we subsequently evaluated the genetic landscape of 11 additional pediatric and adolescent/young adulthood meningioma patients within the Children's Brain Tumor Tissue Consortium. A second individual harboring a YAP1-FAM118B gene fusion was identified within this database. Transcriptomic profiling suggested that YAP1-fusion meningiomas are biologically distinct from NF2-driven meningiomas. Similar to other meningiomas, however, YAP1-fusion meningiomas demonstrated overexpression of EGFR and MET. DNA methylation profiling further distinguished YAP1-fusion meningiomas from those observed in ependymomas. In summary, we expand the genetic spectrum of somatic alteration associated with NF2 wild-type meningioma to include the YAP1-FAM118B fusion and provide support for aberrant signaling pathways potentially targetable by therapeutic intervention.

Disease-associated mosaic variation in clinical exome sequencing: a two-year pediatric tertiary care experience.

Exome sequencing (ES) has become an important tool in pediatric genomic medicine, improving identification of disease-associated variation due to assay breadth. Depth is also afforded by ES, enabling detection of lower-frequency mosaic variation compared to Sanger sequencing in the studied tissue, thus enhancing diagnostic yield. Within a pediatric tertiary-care hospital, we report two years of clinical ES data from probands evaluated for genetic disease to assess diagnostic yield, characteristics of causal variants, and prevalence of mosaicism among disease-causing variants. Exome-derived, phenotype-driven variant data from 357 probands was analyzed concurrent with parental ES data, when available. Blood was the source of nucleic acid. Sequence read alignments were manually reviewed for all assessed variants. Sanger sequencing was used for suspected de novo or mosaic variation. Clinical provider notes were reviewed to determine concordance between laboratory-reported data and the ordering provider's interpretation of variant-associated disease causality. Laboratory-derived diagnostic yield and provider-substantiated diagnoses had 91.4% concordance. The cohort returned 117 provider-substantiated diagnoses among 115 probands for a diagnostic yield of 32.2%. De novo variants represented 64.9% of disease-associated variation within trio analyses. Among the 115 probands, five harbored disease-associated somatic mosaic variation. Two additional probands were observed to inherit a disease-associated variant from an unaffected mosaic parent. Among inheritance patterns, de novo variation was the most frequent disease etiology. Somatic mosaicism is increasingly recognized as a significant contributor to genetic disease, particularly with increased sequence depth attainable from ES. This report highlights the potential and importance of detecting mosaicism in ES.

Samovar: Single-Sample Mosaic Single-Nucleotide Variant Calling with Linked Reads.

Linked-read sequencing enables greatly improves haplotype assembly over standard paired-end analysis. The detection of mosaic single-nucleotide variants benefits from haplotype assembly when the model is informed by the mapping between constituent reads and linked reads. Samovar evaluates haplotype-discordant reads identified through linked-read sequencing, thus enabling phasing and mosaic variant detection across the entire genome. Samovar trains a random forest model to score candidate sites using a dataset that considers read quality, phasing, and linked-read characteristics. Samovar calls mosaic single-nucleotide variants (SNVs) within a single sample with accuracy comparable with what previously required trios or matched tumor/normal pairs and outperforms single-sample mosaic variant callers at minor allele frequency 5%-50% with at least 30X coverage. Samovar finds somatic variants in both tumor and normal whole-genome sequencing from 13 pediatric cancer cases that can be corroborated with high recall with whole exome sequencing. Samovar is available open-source at https://github.com/cdarby/samovar under the MIT license.

Expansion of B4GALT7 linkeropathy phenotype to include perinatal lethal skeletal dysplasia.

Proteoglycans have a core polypeptide connected to glycosaminoglycans (GAGs) via a common tetrasaccharide linker region. Defects in enzymes that synthesize the linker result in a group of autosomal recessive conditions called "linkeropathies". Disease manifests with skeletal and connective tissue features, including short stature, hyperextensible skin, and joint hypermobility. We report a family with three affected pregnancies showing short limbs, cystic hygroma, and perinatal death. Two spontaneously aborted; one survived 1 day after term delivery, and had short limbs, bell-shaped thorax, 11 ribs, absent thumbs, and cleft palate. Exome sequencing of the proband and one affected fetus identified compound heterozygous missense variants, NM_007255.3: c.808C>T (p.(Arg270Cys)) and NM_007255.3: c.398A>G (p.(Gln133Arg)), in B4GALT7, a gene required for GAG linker biosynthesis. Homozygosity for p.(Arg270Cys), associated with partial loss of B4GALT7 function, causes Larsen of Reunion Island syndrome (LRS), however no previous studies have linked p.(Gln133Arg) to disease. The p.(Gln133Arg) and p.(Arg270Cys) variants were transfected into CHO pgsB-618 cells. High protein expression of p.(Gln133Arg) was found, with mislocalization, compared to p.(Arg270Cys) that had a normal Golgi-like pattern. The p.(Gln133Arg) had almost no enzyme activity and little production of heparan sulfate GAGs, while p.(Arg270Cys) only had 17% of wild-type activity. These findings expand the phenotype of B4GALT7-related linkeropathies to include lethal skeletal dysplasia due to more severe loss of function.